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Abstract: We consider a particular case of the two dimensional Blume- 

■ 

Emery-Griffiths model to study the finite-size scaling for a field driven 
• i-H , 

first-order phase transition with two coexisting phases not related by a 

c3 ' 

symmetry. For low temperatures we verify the asymptotic (large volume) 
predictions of the rigorous theory of Borgs and Kotecky, including the 
predictions concerning the so-called equal-weight versus equal-height con- 
troversy. Near the critical temperature we show that all data fit onto a 
unique curve, even when the correlation length £ becomes comparable to 
or larger then the size of the system, provided the linear dimension L of 
the system is rescaled by £. 
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1. Introduction 

First-order phase transitions play an important role in physics. Exam- 
ples range from well known transitions in solid state physics (like ordinary 
melting) to phase transitions in the early universe. For many models de- 
scribing such a transition, exact solutions are not available, and systematic 
expansions such as "weak coupling" or "high temperature series" are of 
limited value. In this situation, the numerical simulation of physically in- 
teresting models has gained more and more popularity. With the increasing 
speed of computers, the statistical errors of Monte Carlo simulations have 
been drastically reduced, but a systematic limitation of this method is the 
finite size of the simulated system. A clear understanding of finite-size 
effects is therefore more and more important in order to interpret the re- 
sulting numerical data. 

The theory of finite-size scaling (FSS) for second-order phase transi- 
tions has a rich and longstanding literature, starting with the pioneering 
work of Fisher [1] and Fisher and Barber [2]. In the last ten years, the the- 
ory of FSS has been extended to first-order transitions. For systems with 
two coexisting phases, three main approaches can be distinguished in the 
literature: the phenomenological renormalisation group approach to FSS 
of Berker/Fisher [3] and Blote/Nightingale [4], the thermodynamic fluc- 
tuation theory for the order parameter distribution iV(s) of Binder and 
co-workers [5,6,7] and the transfer matrix method of Privman and Fisher 
[8], see Refs. [9] and [10] for reviews. 

Recently, the theory of FSS near first order transitions has been put 
on a rigorous basis [11,12] within the context of the Pirogov-Sinai theory 
[13] which applies whenever it is possible to describe the configurations 
of the system in terms of energetically unfavourable contours. Typical 
models to which these methods can be applied are field-driven transitions 
in low temperature spin systems with discrete spins, temperature-driven 
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transitions for Potts models, lattice field theories with double well potential, 
Higgs/confinement transitions in large N lattice gauge theories, and certain 
2d continuum field theories, such as P{4>)2 held theories in the broken phase 
and supersymmetric Wess-Zumino models. 

In the context of the field-driven transitions with two coexisting phases 
considered in this paper, the main result of Ref. [11] can be summarized by 
a formula of the form 

Z(V,h) = L-PMPMv + e -(3f-(f3,h)v^ (i + (y e - L / Lo )) (1-1) 

for the partition function in a cubic volume V = L d with periodic boundary 
conditions [14]. Here L is a constant of the order of the infinite volume 
correlation length, and / m (/5, h) is some sort of metastable free energy for 
the phase m. It is equal to the free energy /(/?, h) if m is stable, strictly 
larger then f([3, h) if m is unstable, and may be chosen as a smooth function 
of h. The finite size scaling of the partition function and its derivatives is 
obtained by expanding / m (/3, h) around the transition point h t . 

In a recent paper [16], we have analysed a controversy stemming from 
the work of Binder and co-workers: within the phenomenological fluctua- 
tion theory for the order parameter distribution Py(s), different formulae 
for the FSS of the order parameter are obtained assuming different ansatze 
[6,7] for the relative normalisation of the two Gaussians. If the normalisa- 
tion is chosen in such a way that both Gaussians have equal weight at the 
infinite volume transition point ht, the resulting formula for the magneti- 
sation in a cube with periodic boundary conditions leads to the prediction 
[6] that the shift of the susceptibility maximum with respect to h t is pro- 
portional to L~ 2d , while a normalisation with equal peak height at ht leads 
to the prediction [7] that this shift is proportional to L~ d . 

To some extent, this controversy had already been resolved in Ref. [11], 
where it had been shown that the predictions of the equal weight prescrip- 
tion agree with those following from equation (1.1) above. Unfortunately, 



equation (1.1) is only proven for very low temperatures (for Ising type sys- 
tems in d = 2, this restricts the temperature to about half T c ), while most 
numerical simulations are performed near T c . It therefore seemed desirable 
to test the predictions following from (1.1) in an actual numerical simula- 
tion in order to definitely settle the controversy between [6] and [7] . 

This was done in [16]. We considered a spin model with spins s x G 
{—1,0, 1} and Hamiltonian 

H = \ \ a M -<?( s y)\ 2 ~ h J2 Sx ' ( L2a ) 

<xy> x 

where the first sum goes over the set of dL d nearest neighbours in V and 

'(«■)={+! E^+i. (1 - 2b » 

It is defined in such a way that boundaries between —1 and or +1 are 
energetically unfavourable, while boundaries between and +1 do not cost 
energy. As a consequence, the low temperature phases of the model consist 
of a phase "minus" corresponding to spin configurations where most spins 
s x are equal to —1 with small islands of or +1, and a phase "plus" 
corresponding to spin configurations which are a mixture of and +1 with 
small islands of —1. The phase minus is stable as long as h is smaller then 
a certain critical value h t = h t ((3), while the plus phase is stable for h > h t . 
Because in the "plus" phase spins can be flipped from to +1 without 
creating any energetically unfavourable boundaries we expect the magnetic 
susceptibility of this phase to be larger than that of the "minus" phase. 
A large susceptibility difference between the phases is desirable because it 
leads to a clear difference between the predictions of the equal-weight and 
equal-height theories of first order finite size effects. 

This model is a special case of the Blume-Emery-Griffiths (BEG) model 
[17], which is the most general model with a three- valued spin and nearest 
neighbour interactions. The general BEG model has a complicated phase 

4 



structure, with up to three coexisting phases. Here we have chosen the 
coupling constants in such a way that the model has two coexisting phases 
at the phase transition point ht, with a large susceptibility difference. 

In [16], we simulated the model (1.2) in d = 2 and determined the shift 
of the susceptibility maximum with respect to h t . For (3 = 0.5, and lattices 
with linear dimensions L from 4 to 30, we found overwhelming evidence for 
the equal weight prescription, as to be expected from the rigorous results 
of Ref. [11]. We summarize the main results of Ref. [16], as well as certain 
extensions (in particular to (3 = 0.45 and (3 = 0.47) and refinements in 
Section 3. 

In [16], we did not analyse, however, other finite-size scaling predic- 
tions following from equation (1.1), such as the volume dependence of the 
susceptibility-maximum, or the volume dependence of the susceptibility at 
the transition point h t . In view of the results of Billoire and co-workers 
[18], who found that the corresponding asymptotic scaling form for two di- 
mensional Potts models was only reached for systems of linear dimensions 
L as big as 5 to 10 times the correlation length £ of the infinite system, it 
seems highly desirable to check more finite-size scaling predictions for the 
model considered above. This will be done in Section 4. 

As we will see, the asymptotic FSS-form obtained from equation (1.1) 
will again only be reached for fairly large volumes. While these volumes are 
easily accessible for (3 = 0.5 and (3 = 0.47 (the corresponding correlation 
length £ is 2.19238... and 4.35091...), it is not accessible for temperatures 
near T c (we did do a third simulation for (3 = .45, corresponding to £ = 
13.5096...). For these temperatures, however, a scaling ansatz using the 
length scale £ to rescale all lengths in consideration, should be reasonable. 
In fact, we will see in Section 5, that this ansatz allows to fit all data onto a 
unique curve which describes the crossover between first and second order 
phase transitions. 
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2. The Model and its Relation to the standard Ising Model 

In the next three sections, we will analyse the FSS for the model (1.2). 
We recall that the low temperature phases of the model consist of a phase 
"minus" corresponding to spin configurations where most spins s x are equal 
to —1 with small islands of or +1, and a phase "plus" corresponding to 
spin configurations which are a mixture of and +1 with small islands of 
— 1. As we will see below, it can be shown that for f3 > f3 c , where (3 C is the 
critical temperature of the Ising model, the phase minus is stable as long as 
h is smaller then h t = — /3 _1 arsinh(l/2), while the plus phase is stable for 
h > ht. At h = ht, the infinite volume magnetisation jumps from a value 
m_ = m-(P) to a value m + = m + ((3) > m_. 

In order to compare the predictions of FSS theory with numerical data, 
one needs the numerical values of the infinite volume magnetisation m± and 
susceptibility x± & t the point h = h t ±0. As shown in [16], these constants 
can be obtained as follows: Let Z(V,P, h) be the partition function of the 
model (1.2), and let Zj(V,/3,h) be the partition function of the standard 
Ising model. Then Z and Zi are related by the formula 

Z(V, p, h) = e W,h)-h)pv Zi ^ ^ ^ (2 _ la) 

with 

(jl(J3, h) = (1/2/3) He^ + e 2f3h ) . (2.1b) 

As a consequence, we obtain the free energy /(/?, h) of the model (1.2) 
in terms of the free energy //(/?, n) of an Ising model with magnetic field fx 

f(p, h) = f T (p, fj,(p, h)) + h- »(p, h) . (2.1c) 

In a similar way, it is possible to relate the correlation length £ and the 
interface tension a of the model (1.2) at h = h t to those of the Ising model 
at zero field. One finds 

m = ZiiP) and a(P) = aj(p). (2.2) 
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Recalling that phase transitions are associated with singularities of the 
free energy, we can read off the phase diagram of the model (1.2): It has 
the same critical temperature as the Ising model, while ht is just the point 
where /a(0, h) is zero, i.e. 

1 ■ , /, 0.48121183... 
h t = --arsmh(l/2) = . (2.3) 

Differentiating equation (2.1c) with respect to h and taking the limit h — > 
ht ± 0, we finally obtain m± and x±- I* 1 terms of the magnetisation mj = 
mi((3) and the susceptibility xi = Xi(P) °f the Ising model at [i = +0 one 
gets 

m± = ~e ph * ±^(2-e^*)m 7 (2.4a) 
X± = ^(/3±/3m 7 )e 3 ^ + 1(1 + e 2 ^)\j . (2.4b) 
Note that the difference x+ ~ X- does not involve the constant xii 

X+~X- = (Pmi)e spht • (2.4c) 

We have simulated the model (1.2) in d = 2, where m/(/3) and £/ = 
1/2(7/ are exactly known [19], 

mi (j3) = (1 - (sinh2/3)- 4 ) 1/8 (2.5a) 
^ (/3) = (4/3 + 21n(tanh /3 ) ) " 1 for > f3 c 

= -(2/3 + ln^anh/?))" 1 for (3 < C . (2.5b) 

In order to compare FSS theory to numerical data, we will also need x+ an d 
X-- Unfortunately, the susceptibility Xi(P) °f the Ising model is not exactly 
known. We therefore used the low-temperature series for Xi(P) found in 
[20] together with the method of Pade approximants [21] to calculate X/(/3). 
As in [22], we embodied the knowledge of the exact transition point and 
critical exponent in our approximation. The result is 

Xi{0)=A0u* |l + 5>n« n j ^(l-6u + u 2 )^l + J^q n u n Jj 

(2.6) 



where u = e _4/3 , and p n , q n are given in table 1. The factor 1 — 6u + u 2 
has a simple zero at (3 = f3 c , ensuring that xi diverges as (f3 — /3 C ) 7 / 4 . 
It is amusing to note that a Pade approximant of the form (2.6) delivers 
the exact formula (2.5a) for m; if used on the strong coupling series for 
magnetisation. 



n 


Pn 


q n 


n 


Pn 


q n 


1 


-11.422 


-9.9941 


4 


-0.83867 


-41.824 


2 


42.042 


27.744 


5 


-4.9026 




3 


-49.064 


-5.1210 









TABLE 1: Low temperature Pade coefficients for the susceptibility of the 
Ising model. 



8 



3. Equal Weight versus Equal Height 

In order to obtain the FSS of the order parameter distribution -FV(s) 
of an order parameter s (e.g. s = y Yl s x is the magnetisation per lattice 
site for spin models ) in a cubic system with periodic boundary conditions 
and volume V = L d , Binder and Landau approximated the probability 
distribution Pv(s) by a sum of two Gaussians, centered around the infinite 
volume magnetisation m_ and m+. As for the relative normalisation, in 
the original version [6] they added up two normalised Gaussians so that at 
the transition point ht each of the two coexisting phases contributes with 
the same weight. Following this "equal weight" -rule [6] one gets 



Pv{s) * 1 J_ e -^^— + _J_ e -^^x^- e Ps(h-h t )V 

(3.1a) 

where x± are the infinite volume susceptibilities at the point ht — and 
h t + 0, respectively, and iV = N((3, h) is chosen in such a way that Pv(s) 
is normalised. Led by mean-field arguments they revised this version for 
the relative normalisation, adding up the two Gaussians in such a way that 
the peak heights are equal at the transition point ht- Following this "equal 
height"-rule [7] Pv(s) turns out to be 

JV(a) = ^ e~ VP ^— + e~ VP ^^- e^ h - h ^ v . (3.1b) 

It is an easy exercise to calculate the FSS for the finite volume magnetisa- 
tion m{V,h) = J sPv(s)ds following from (3.1a) and (3.1b). For the point 
h x (V) where the susceptibility x(V, h) = dm(V, h)/dh has its maximum, a 
lengthy but straight-forward calculation gives the FSS predictions 

and 
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h x (V) = h t 



Hx+/x-) i 

2/3(m+ -m_) V 




X+ ~ X- 1 
f3 2 (m+ — m_) 3 



+ 0(V" 3 ) , (3.2b) 



respectively. We emphasize that the shift predicted by the "equal weight" - 
rule [6] is proportional to 1/V 2 while the shift predicted by the "equal 
height" -rule [7] is proportional to 1/V. (For symmetric phase transitions 
both rules are equivalent and predict that h x (V) = h t .) 



Fig. la: The position h x (V) of the susceptibility maximum for (3 = 0.45. 
Open points are Monte Carlo results, solid points are exact results ( see ap- 
pendix). The solid line is the theoretical prediction (3.2a), while the dashed 
line is the prediction (3.2b). 

We have tested the predictions (3.2) for three different temperatures, 
namely for (3 = 0.45, (3 = 0.47 and (3 = 0.5. In our simulation we first 
generated the probability distribution -Py(s) for L = 4, 5, 12, 14, 18, 30 
at h = h t using a multimagnetical version [23] of the heat bath algorithm. 
Applying the reweighting technique of Ferrenberg and Swendsen [24], we 
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obtain the probability distribution -Py(s) and hence the magnetisation, 
m(V,h) = J2 s Pv(s) and the susceptibility, h) = PV{J2 s 2 Pv(s) — 
sPv(s)) 2 }, for h in the vicinity of ht- In order to use the vector structure 
of the CRAY X-MP and Y-MP, we calculated 64 independent copies of the 
system in parallel, with 15 million sweeps for each of them. Grouping 8 
of them together to save memory, we obtained 8 statistically independent 
data sets for each L, and therefore 8 statistically independent estimates 
for h x (V). Statistical error estimates are obtained from these data in the 
standard way. 



Fig. lb: The position h x (V) of the susceptibility maximum for (3 = 0.47. 

We have plotted the theoretical curves (3.2a) and (3.2b) together with 
our numerical data in Figs. la-c. We recall that the coefficient in (3.2a) is 
known exactly, while the coefficients in (3.2b) have been obtained with the 
help of a Pade improved low-temperature expansion. Within the resolu- 
tion of Fig. 1, the errors for the numerically determined points h x (V) are 
invisible, except for f3 = 0.45. Fig. 1 shows overwhelming evidence for the 
equal weight prescription. 
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Fig. lc: The position h x (V) of the susceptibility maximum for (3 = 0.50. 

Convinced of the "equal weight" -rule we tried the new method to de- 
termine the transition point ht of first-order phase transitions recently pro- 
posed by C. Borgs and W. Janke [25]. This method just uses the fact that 
both phases should have equal weight at h = ht- To be precise, one defines 
a finite volume transition point hEw(V) as the point where 

Pv{s)= p y^- ( 3 - 3 ) 

Here s m i n , s_ < s m i n < s+, is the point s where the distribution Pv(s) 
has its minimum. According to Ref. [25], it is expected that the system- 
atic errors \h,Ew(V) ~ h t \ are exponentially small. Here they are so small 
that hEw{V) agrees with h t within the statistical errors as soon as L/£ is 
noticeably bigger then one (see Table 2 at the end of this paper). 
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4. Asymptotic First-Order FSS-Scaling 

In this section we test the FSS predictions following from (1.1) by 
Taylor expanding f±(/3, h) around h = h t , 

f±(P,h)=f(P,h t )-m ± (h-h t )-^(h-h t ) 2 + ^(h-h t f + --- . (4.1) 

For the maximum x max (^) of the susceptibility, its position h x (V) and the 
magnetisation m(h x (V)) at the point h x (V), the FSS predictions following 
from (1.1) are 

h ^ = h ' + fZl-™-)^ + °^ (4 ' 2 » 

(u m\\ rn+ + m- , 3 X- ~ X+ 1 , n ( 1 \ < A n\ 



Xma X (^)=/3l I V + +°{y) 



(4.4) 



It is also interesting to analyse the FSS of the susceptibility and magneti- 
sation at the infinite volume transition point ht (for systems where ht is 
not known exactly, one could use the point Iiew(V) introduced in [25]). 
Due to the exponentially small error bound in (1.1), the corrections to the 
known leading order FSS predictions are exponentially small 

m_|_ + m_ 



X(h) 



m(h t ) = m ++ m ~ (i + (Ve- L / L °)) (4.5) 
P ( m+ ' 2 m ~ ) 2 V + X++ 2 X ~ ] (l + 0(Ve- L ^)) . (4.6) 



Note that the error bounds in (4.5) and (4.6) are only bounds. For the model 
considered here, the value of m(h t ), e.g., is in fact exactly (m + — m_)/2, 
due to (2.1) and the fact that the finite volume magnetisation with periodic 
boundary conditions and no external field is zero for the standard Ising 
model. The corrections for x(/it), on the other hand, are expected to be 
asymptotically proportional to L a e~ L ^, with a constant a not necessary 
equal to the worst case bound (4.6). 
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In addition to the FSS of the susceptibility and magnetisation, we also 
analyse the FSS of the cumulant Uy introduced by Binder [5]. It is defined 

as 

u (M 4 ) c _ ((M-(M)) 4 )-3((M-(M)) 2 ) 2 

V 3{M 2 ) 2 c 3((M- (M)) 2 ) 2 

Here (-) c denotes connected expectation values and M = J2 x ev Sx - By 

the inequality (F 2 ) > (F) 2 (applied to F = (M - (M)) 2 ) the cumulant 

Uy cannot exceed the value 2/3. In the high temperature region, where 

the asymptotic probability distribution of M is gaussian, Ul goes to in 

the infinite volume limit for all h, while in the low temperature region 

considered here 



Uv(h t ) = I 



1- 4 (* + + X-) 1 + Q(J-)' 

P(m+-m.) 2 V W 2> 



(4.8) 



implying that Uy goes to 2/3 at the first order transition point h = ht if 
V — > oo. At the second order transition point h = h t , j3 = j3 c the cumulant 
Uy is expected to go to a non-zero limit strictly smaller then 2/3. 

In most theories ht is not known. Instead of Uy(ht), it is therefore 
more natural to consider the maximum Uy &yi of Uy as a function of h, and 
its position hu{V). Again, the FSS of these two quantities can be obtained 
from (1.1). One obtains 

^ y »-"- + ^ + + :i + ')3 ^ +0( ^ ) (49) 



v 3 



1- 4 (*+ + *-) 1 + Q(J-) 
(3(m+-m-) 2 V W 2) 



(4.10) 



Even though (1.1) is rigorous for large /?, and probably correct for 
all (3 > P c provided V is large enough, it is not clear where the large V 
asymptotics for the first-order FSS predictions (4.2) through (4.10) sets in. 
In fact, warned by the results of Billoire and co-workers [18], we expect 
that the asymptotic FSS sets only in once where 3> might mean 

five to ten times larger. 
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We have tested the above finite-size scaling predictions for three dif- 
ferent temperatures, namely for (3 = 0.45, (3 = 0.47 and (3 = 0.5 (cor- 
responding to f = 13.5096..., f = 4.35091... and f = 2.19238...). In our 
simulation we measured the quantities h x (V), m(h x (V)), Xmax(^), x(ht), 
hu(V), and U v (h t ) for L = 4, 5, 12, 14, 18, 30, see Section 3 above 

for the details of the simulation. The results are given in Table 2a-c. 

We show both the theoretical prediction according to (4.2) through 
(4.10) and our numerical results for h x (V), hu(V), m(h x (V)), x(ht) and 
Uv(ht) in Figs. 2-5 below. Monte Carlo data are represented by open 
circles and open squares, and theoretical predictions are represented by full 
or dashed lines. 

As a first point we want to stress that our theoretical calculations 
predict that both h x (V) — h t and hu{V) — h t go to zero like 1/V 2 , a 
peculiar feature of an asymmetric first order transition with two coexisting 
phases 1 . Indeed, our numerical data conform with this prediction, see 
Fig. 2a. Since the last points in Fig. 2a are barely distinguishable, we 
show a blow up of this figure in Fig. 2b. Unfortunately, the statistical 
errors for the corresponding data points are so big that our data for h x (V) 
and hjjiV) agree for the last five lattice sizes L =11, 30, see Table 2. 
We therefore decided to present data points from an exact transfer matrix 
calculation in Fig. 2b. As can be seen from the figure, these data are in 
perfect agreement with the theoretical predictions (4.2) and (4.9). 

Next, we point out that the volumes for which the asymptotic FSS 
according to (4.2) through (4.10) sets in depends strongly on the considered 
quantity. While the asymptotic regime corresponds to L/£ > 1 — 2 for 

1 More generally, h x {V) - h t ~ 1/V 2 and h v (V) - h t ~ 1/V 2 if the 
number N- of phases stable for h < h t and the number N + of phases 
stable for h > h t are equal. If ^ N_ the shift in h is proportional to 
1/V. 
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Fig. 2a: The finite-size transition points h x (V) (dashed line) and hjj(V) 
(solid line). 
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Fig. 2b: h x (V) and h v (V) (blow-up). 



17 



Fig. 3: The magnetisation at the susceptibility peak 
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Fig. 4-' The susceptibility x(ht) an d its asymptotic behaviour. 
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Fig. 5: The cumulant Uy(h t ) and its asymptotic behaviour. 
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he(V) and h v (V), and to L/£ > 2 - 3 for m(h x (V)), both x(h t ) and U v (h t ) 
require much larger volume (L/£ > 4 — 5 and L/£ > 6 — 8, respectively). 
This conforms to the results of Gupta and Irback, who observed similar 
effects for a symmetric field driven transition [26] . 

A last point of interest is related to the determination of the transition 
point ht as discussed at the end of the last section [25,27]. In addition to 
the point hEwiY) defined in (3.3), one may consider "single peak" mag- 
netisations and susceptibilities m±(V) and x±(V) defined as 

and 

X±(V)=(3v(J2s 2 P$(s)-(J2 sP v( s )) 2 ) > ( 4 - 12 ) 

where Py is the normalised probability distribution of the left and right 
peak, respectively. More explicitly, 

m .(v)= sp y^/ E p y^ 

m+(V)= sP v( s )/ E p v( s ) 



(4.13) 



and similarly for x±(V)- 

The corresponding data can be found in Table 2a-c as well. An in- 
spection of these data for (3 = 0.50 shows an interesting feature. For small 
system sizes, x±(Y) are smaller then the infinite volume values x±i an< i 
m+(y) > m + and while m_(V) < m_. As the system becomes larger, all 
these quantities approach the infinite volume value, and in fact overshoot, 
before they finally approach the infinite volume limit from the other side, 
see Fig. 6. We have observed the same effect for even lower temperatures 
(where our exact transfer matrix calculation can easily reach the asymp- 
totic regime), and we expect that it occurs for all temperatures (3 < (3 C 
provided the lattice is chosen large enough. 
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We think that this effect is due to the competing influence of two 
finite-size effects for the probability distribution Py(s). One of these ef- 
fects concerns the fact that the single peaks Py obtained by just cutting 
the original distribution at the point s m i n contain events which actually 
should be considered part of the other phase. A minute of reflection shows 
that this effect decreases the susceptibility, and shifts the magnetisations 
m±(V) in the direction observed for small lattice systems (if one assumes 
that the original distribution is of the form (3.1a), this can be shown rigor- 
ously). The second, competing effect comes from tunneling events. These 
events tend to shift the magnetisations towards each other, and to increase 
the susceptibility. Since overlap effects between the two peaks die out ex- 
ponentially with the volume V = L d , while tunneling effects die out expo- 
nentially with L d_1 , the competition of these two effects seems to naturally 
explain the observed overshooting for m±{V) and x±(Y)- 



Fig. 6: The volume dependence of the "single peak" magnetisation m + (V) 
for (3 = 0.5. The dashed line shows the infinite volume value. 
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5. The Crossover from First- to Second-Order FSS 

We have seen that the formulae of sections 3 and 4 break down when 
the correlation length £ becomes comparable with the lattice size L. In this 
section we try to explain finite size behaviour in this large £ region. 

In Fig. 7 we show the susceptibility at ht as a function of lattice size for 
a range of inverse temperatures f3 in the neighbourhood of (3 C . The results 
for the smaller lattices (up to size 14 x 14) are taken from exact calculations 
of the magnetisation distribution (see appendix). The susceptibilities on the 
larger lattices are Monte Carlo results, either from our own simulations with 
the three-state model (1.2) or derived (via the formulae in section 2) from 
simulations of the two dimensional Ising model [23,29]. What we see is that 
on small lattices the curves are all rather similar, independent of whether 
(3 is above, below or at the critical point. This is indeed reasonable, on a 
finite lattice all physical quantities are analytic functions of the coupling 
constants so there can be no sharp distinction between the first order phase 
transition above (3 C , the second order phase transition at j3 c or the rapid 
crossover (without a phase transition) seen below (3 C . 

This suggests that we should try to understand this region by using 
the results of the finite size scaling scaling theory in the neighbourhood of 
a second order phase transition. (For reviews see [10,28]). We expect this 
theory to apply whenever the correlation length £ and lattice size L are 
both large compared to the lattice spacing. 

The basic hypothesis is that when the correlation length is large enough 
thermodynamic properties are no longer dependent on the lattice spacing 
itself. By this is meant that instead of being a function of the three param- 
eters f3, h and L the free energy and correlation length become functions 
of only two variables, i.e. 
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Fig. 7: x(ht)/V as a function of lattice size. Solid points are exact results, 
open points are Monte Carlo results. 
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f ing (/3, h, L) « L- d /(tL 1 ^, HL^+M") (5.1a) 
and 

fc, L) « L£(tL 1/u , HLW V ) (5.1b) 

where the thermodynamic variables are t = (/3 C — (3)/ (3 and H = h — h t , 
and 7,/5 and z/ are the critical exponents of the model (1.2), which (as 
explained in section 2) are the same as those of the Ising model, namely 
7 = 7/4, (3 = 1/8 and v = 1. (It is unfortunate that the traditional symbol 
for inverse temperature and that for the critical exponent are both /3, we 
hope that the context will always make it clear in which sense (3 is to be 
understood.) 

Differentiating the free energy leads to a scaling form for the magnetic 
susceptibility 

X(P, h, L) » L^xitr* 1 '", HlSi+M"). (5.2) 

(The "hyperscaling" relationship 2/3 + 7 = dv has been used to simplify 
this expression.) A simple test of the scaling form (5.2) comes from our 
data at h = h t . The finite size scaling prediction is 

X (/3,^,L)/L^^x(^,0), (5.3) 

i.e. that all data should lie on a universal curve 

if x(/3, h t ,L)/L^/ v is plotted 
against tL x l v '. The results are shown in Fig. 8. In place of tL 1 l v we have 
used the somewhat more physical coordinate 

(4/3 + 21n(tanh/3))L 1 / iy = (L/£(J3, h t , oo)) 1/v for [3 > (3 C 

= -2(L/£(f3,h u oo)) 1 /» for (3 < (3 C .(5.4) 

This is possible in this model since we have an exact expression (2.5b) for 
the correlation length. (The factor —2 in the symmetric phase is there to 
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ensure that the ^-coordinate is an analytic function of (3.) We see that 
the data do indeed fall on a universal curve. (Some deviation is seen at 
(3 = 0.50, which is to be expected since the correlation length is only 
& 2.19). The curves predicted by the ansatz (3.1a) for first order phase 
transitions are also plotted in this figure. We see that the predictions are 
reliable as long as L/£ > 4. The asymptotic results calculated in section 3 
allow us to say more than the scaling hypothesis (5.1), because they give 
us the large x asymptote of the scaling function x(x, 0). 



Fig. 8: A scaling plot of x(h t ). Exact results are shown as filled symbols, 
Monte Carlo results as open. Only data for L > 3 are plotted. The curves 
show the expected large volume behaviour. The dashed curves are the large 
volume asymptotes for (3 = 0.30,0.47 and 0.50. The solid curves are the 
large volume asymptotes at (3 = (3 C ± e. 

The scaling form (5.2) also has implications for h x , the position of the 
susceptibility maximum. 

h x (0, L) - h t {(3) w L-^ + M v H(tL 1 / v ) (5.5) 

Taking higher derivatives of the free energy leads to scaling expressions 
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for any cumulant of the magnetisation distribution, in particular to a scaling 
expression for the Binder parameter Uy 

U v (p, h, L) w \J{tL x l\ HL<-i + M v ) (5.6a) 

and its maximal value 

U^ ax « £/ max (£L 1 / I/ ). (5.6b) 

The magnetic field at which Uy reaches its maximum is given by a rela- 
tionship of the form (5.5). 

The scaling predicted in (5.6b) is tested in Fig. 9. As in section 4 we 
find that larger lattices are needed to see finite size scaling for Uy than 
for x- Scaling only works well on lattices with L larger than about 5, on 
smaller lattices the value of £/"™ ax lies noticeably below the scaling curve. 



Fig. 9: A scaling plot ofUy a *. Exact results are shown as filled symbols, 
Monte Carlo results as open. Only data for L > 5 are plotted. The curve 
is the large volume asymptote (4-10) at (3 = (3 C + e 

In conclusion we see that there is a considerable region of overlap be- 
tween the domain of validity of the first order finite size scaling relationships 
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presented in earlier sections and the scaling region results presented here. 
The first order results require that L be large with respect to the correla- 
tion length, but make no demands on what the correlation length must be, 
while the results of this section require that L and £ are both large, but 
do not restrict the ratio L/£. The domain of mutual validity is therefore 
L >> £ >> 1, while in the entire range L > 1 we can always apply at least 
one of the formulae. 
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Appendix - Exact Calculations 

Naively one might think that exact calculations of the partition func- 
tion for our model would be completely impractical, since the number of 
configurations that need to be considered is 3 L . However by carefully 
dividing the lattice up into smaller sub-lattices it is possible to construct 
algorithms in which the amount of calculation grows exponentially with 
L d_1 rather than L d , which gives spectacular savings in time, especially 
when d = 2. 

The essential insight is that the partition function for a lattice spin 
model can be calculated recursively, see [30,31]. 

Suppose that the partition function Zn for the N site sub-lattice shown 
at the top of Fig. 10 is a known function of the spins on the boundary, (a\ 
to oil)- The spins a in the interior of the lattice (the gray region) have 
been summed over. Periodic boundary conditions in the x direction are to 
be understood. 

Zn(&1, ■ ■ ■ ) &L', &L+1, • • 4 , Cj-l! <7j, • • • ) &2l) 

= y%xp(-/?.ffjv(or,cri,- • -,(T2l)) (A.l) 
a 

The lattice size is next increased by one by adding the new spin <jj. The 
energy difference, 5H, between the new lattice and the old lattice depends 
only on the boundary spins, not on any of the interior spins (in the example 
shown the energy change depends only on the three spins <jj and u'-, 
and it is also true for the special cases of the first and last spin in a row 
that only boundary spins are involved in the energy change). 

Finally we can find the partition function for the larger lattice by 
summing over aj. Algebraically 

ZN+l((Tl, ■ ■ ■ , &L', 0~L+li ■ ■ ■ , Cj', + ■ ■ ■ , U2l) 
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Fig. 10: The procedure for adding a spin to the lattice when calculating the 
partition function recursively. (See the text for a detailed explanation.) 
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= J^exp (-(3 5H(a L+lr --,a L ;a' j )) 

xZ N (a 1 , ■ ■ ■ ,a L ; a L+1 , ■ ■ ■ , aj-i, aj, ■ ■ ■ , o^l) (^- 2 ) 

The recursion can be started from the fact that Zq (the partition func- 
tion of a lattice with no spins) is 1. 

In our implementation of this procedure the amount of calculation 
grows as q 2L and the amount of storage needed as q L for a q-state spin 
model. The basic reason that the calculational burden grows like q 2L is 
that the partial partition function Zjv depends on the 2L boundary spins, 
which means that its arguments can take q 2L different values. The number 
of arguments would be greatly reduced if the boundary conditions in the 
y direction allow us to sum over the spins in the first row (the spins a± to 
0"^). Examples of boundary conditions which permit such a sum are either 
free boundary conditions or a fixed boundary configuration. These would 
reduce the computational time to a value oc q L . Even though they are not 
the best boundary conditions for exact calculations we have used periodic 
boundary conditions in both directions to make comparison with Monte 
Carlo results possible. 

The close relationship of our model to the Ising model means that we 
can calculate partition functions for the Ising model which has q = 2 and 
then transform to the q = 3 model at the end, meaning that computation 
grows as 4 L instead of 9 L . 

Even though the calculation grows exponentially we were able to han- 
dle surprisingly large lattices, with sizes of up to 14 x 14, without needing 
a supercomputer. 

Once we have found the multiplicity of states Ni(E,n-) for the Ising 
model it is a simple matter to use the binomial distribution to calculate 
the corresponding multiplicities for our model (1.2) 

N(E, n_, n ) = N^E, n_)— (A3) 

no ] .{L a — 7i- — no)' 
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Here E is the number of broken bonds and n_, no and n + are the number 
of spins of a given type. We have a simple relationship of this type be- 
cause boundaries between spin and spin +1 entail no energy cost while 
the energy cost for boundaries between spin —1 and either of the other 
two spin types is the same in both cases. It is a simple matter to derive 
relationships between partition functions and magnetisation histograms of 
the two models by summing (A. 3) over (some of) its arguments. 
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TABLE 2a: Monte Carlo results for (3 = 0.45 



35 



L 


hJV) 


mihJV)) 

\ A. \ / / 




x(^)/^ 




A 

4 


i ni snon/ c;n^ 

-i.uiouyy^oy j 


ooi 7r>ftQ/'77 N \ 
-.zyl / Uoo(^ / / J 


1 Q70ftQ1 / K 1 ^ 

.io /yooi^oi J 


.io /oo4t>(^ob J 




r 




i nil a £ i / q /i ^ 
-l.Uzl4ol(o4J 


ongnft7/i /ftn\ 
-.zyoUD(4(Dyj 


i gi 7ft 1 1/'£ft^ 
.lol ( DlZ^OD J 


i g i 7nni / /i r\ 
.lol < Uyi(4b J 




ft 



-l.Uzzo4U^o * j 


Qm /i ftfi q ^n"! 
-.oU14ooo(^ 1 UJ 


.1 ( 1 oZooyoo) 


1 77^m 7^ft1 
.1 I lOKji. ({pi. ) 




7 


-l.UZoloZ^Zo ) 


-.oUot)UU4^ / /J 


.1 1 44000(^00 j 


1 7/1 A 7r>Q z' C^Q^ 

.1 / 44 1 Uoyoo ) 




Q 

o 


-l.UzooUU^Z I ) 


-.oU4ozt)U^ol J 


1 lOO A Oft ^7^ 
.1 / ZZ4ZD^0 / J 


.1 / ZZOO / (^0 1 ) 




9 


-1.023655(16) 


-.3057188(65) 


.1705276(51) 


.1705242(51) 




10 


-1.023683(14) 


-.3063460(64) 


.1692254(71) 


.1692216(67) 




11 


-1.023747(20) 


-.3068320(83) 


.1682018(49) 


.1681995(46) 




12 


-1.023794(23) 


-.3071837(68) 


.1674027(69) 


.1674017(67) 




14 


-1.023813(14) 


-.3076732(72) 


.1662249(70) 


.1662241(71) 




18 


-1.023858(09) 


-.3082108(91) 


.1649093(29) 


.1649093(29) 




30 


-1.023846(06) 


-.3087157(58) 


.1636313(95) 


.1636305(96) 




L 




jymax 


Uv(ht) 










4 


-l.UzUzUo^oy ) 


fti nn/iofn^ 

.DIUU^IZ^UO J 


.ouyo i oy io j 










c; 

D 


-l.UzzoZo^o4 ) 


fti Q7fif|('nzL N l 


fti QftOftYnft"! 
.DiyOZD^UD J 













i noQnft7/Q7 N \ 
-l.UzoUb i [o ( ) 


ftofi^oftf'nft'i 

.DZDoZD^UD ) 


ftOftQ/IOfl 1 
.DZDZ4Z^11 J 










7 

< 


1 (T")Q /I 1 Q fOQ'\ 

-l.Uzo41o^zo J 


.Doizuy i^uo j 


ft^i i fti ('nft'i 

.OOI 1D1 yUD 1 










g 
O 


-L.UZovooyZ i ) 


ft^i 9/1 f 

.OODIZ^UD J 


ft^^i mi'ns'i 

.DOOlUO^Uo J 










9 


-1.023740(16) 


.638356(05) 


.638348(06) 










10 


-1.023738(14) 


.641133(07) 


.641119(06) 










11 


-1.023784(20) 


.643530(05) 


.643523(05) 










12 


-1.023820(23) 


.645639(06) 


.645636(07) 










14 


-1.023828(14) 


.649127(06) 


.649124(08) 










18 


-1.023863(09) 


.654133(03) 


.654133(04) 










30 


-1.023847(06) 


.661167(06) 


.661162(09) 










L 




m+(V) 


m_(V) 


X+(V) 


X- 


-(V) 


4 


-l.UZ4UoZ^oU J 


.303542(24) 


-.922636(09) 




.2300(01) 


.1273(01) 





i r>OQOftn/'Qc;\ 
-LUzoyby^oO J 


.296404(16) 


-.915421(15) 




.2807(01) 


.1770(01) 


ft 



i noQ g i ft z' /i c\\ 
-l.UzoolD(4y J 


.291949(99) 


-.909913(99) 




.3280(11) 


.2320(10) 


7 


-l.UZoozZ^ZZ J 


.288400(97) 


-.906357(92) 




.3783(12) 


.2830(12) 


g 
O 


-l.Uzooo I [ZD ) 


.285761(73) 


-.903962(76) 




.4281(11) 


.3297(12) 


9 


-1 023896(17) 


.283968(58) 


-.902061(63) 




.4727(12) 


.3761(13) 


10 


-1.023843(15) 


.282669(71) 


-.900775(64) 




.5146(18) 


.4171(17) 


11 


-1.023858(19) 


.281679(38) 


-.899877(42) 




.5533(10) 


.4536(12) 


12 


-1.023872(23) 


.281054(34) 


-.899190(31) 




.5859(11) 


.4874(11) 


14 


-1.023857(14) 


.280179(28) 


-.898372(28) 




.6433(13) 


.5400(15) 


18 


-1.023873(09) 


.279668(27) 


-.897689(30) 




.7059(19) 


.6123(16) 


30 


-1.023848(06) 


.279645(23) 


-.897639(17) 




.7451(10) 


.6528(09) 


oo 


-1.02385495.. 


.27967484.. 


-.89770883.. 


.7378 


.6433 



TABLE 2b: Monte Carlo results for (3 = 0.47 
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On/1 K1Q K f r >'7\ 

.zU4oooD^z 1 ) 
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nfti fi/ift/'on^ 
-.yolo4o^zU J 


QnQfti (\i I 


ononc;/io/'Qn > \ 

.zuzyo4z(^ou j 


ononQso ( 'i a \ 
.ZUzyooz(^o4 J 




Q 

o 


-.\)uZlZ(j{6Z ) 


on/i ni ft 1 / qq\ 
-.o(j4yLDL{Z6 ) 


oni gftog / /1o^ 
.zL)1odoo(4z J 


oni g £ftQ ( acC\ 
.zU1oodo(4d J 




9 


-.962230(32) 


-.3057789(67) 


.2011162(70) 


.2011113(60) 




10 


-.962293(23) 


-.3064065(53) 


.2005801(45) 


.2005767(48) 




11 


-.962335(20) 


-.3068779(35) 


.2001908(63) 


.2001885(60) 




12 


-.962391(17) 


-.3072287(51) 


.1998887(67) 


.1998882(69) 




14 


-.962375(06) 


-.3076952(53) 


.1994660(22) 


.1994642(24) 




18 


-.962421(14) 


-.3082342(72) 


.1990105(36) 


.1990105(37) 




30 


-.962427(08) 


-.3087290(42) 


.1985327(29) 


.1985326(26) 




L 


h v (V) 


^max 

V 


U v iht) 








4 




ftoonnft/i 
.ozzuyo4^oo J 


ftOI 1 A KQ 1 1 1 
.DZl / 40o(^ll J 








D 


nfti 1 nfi ( aq\ 
-.yt>llUo(^4o J 


ftQOQOI 1 /OC^ 
.DoZoZll^ZO J 


ftQOftQnn/'i o"\ 
.oozooUU^lz J 











nfti 7ft 1 /q /i ^ 
-.ybl f Ol(o4 J 


ftong £Q 1 /q 1 ^ 
.boyoool(ot> J 


.ooy( ioy{0y) 








7 

< 


nftonQ/i / on - ! 
-.yozUo4^zU J 


ft/l/1Q71ft/'1ft N \ 


ft/i/ifiO/in/'ni^ N \ 
.044oz4U(^U0 j 








g 
O 


nftooon/QO^ 
-.yozzoU^oz ) 


ft/I QftOftl /'OQ\ 
.04o0Z0l^ZO J 


ft/i fiftnfti /nn^ 
.04ooUol^Uy J 








9 


-.962298(32) 


.6515315(27) 


.6515180(07) 








10 


-.962337(23) 


.6538180(23) 


.6538081(06) 








11 


-.962365(20) 


.6556541(31) 


.6556476(04) 








12 


-.962413(17) 


.6571408(20) 


.6571405(03) 








14 


-.962387(06) 


.6593655(08) 


.6593586(03) 








18 


-.962425(14) 


.6620330(08) 


.6620329(07) 








30 


-.962428(08) 


.6649091(06) 


.6649074(18) 








L 


h E w(V) 


m+(V) 


m_(V) 


x+(v) 


X-(V) 


4 


nfto k 1 7/ /i r\\ 
-.yozol I (^41) J 


.325905(14) 


-.944591(07) 


.20538(04) 


.09385(03) 





nfto/in/i / xi 
-.yoz4y4(^4o J 


.323234(09) 


-.941739(09) 


.23116(02) 


.11961(03) 


ft 



nftO/i -\-\(oa\ 
-.yoz411^o4J 


.321985(05) 


-.940006(11) 


.25043(05) 


.14281(06) 


7 


nftoogft/on^ 
-.yozooO^zU J 


.321226(27) 


-.939220(26) 


.26679(35) 


.15947(34) 


g 
O 


nfto/i /in/QO^ 
-.yoZ44U^oz J 


.320767(09) 


-.938869(08) 


.27967(08) 


.17092(05) 


9 


- 962428(32) 


.320633(17) 


-.938654(18) 


.28755(27) 


.17991(31) 


10 


-.962423(23) 


.320555(16) 


-.938597(13) 


.29333(28) 


.18540(22) 


11 


-.962423(20) 


.320557(10) 


-.938592(14) 


.29611(12) 


.18939(14) 


12 


-.962454(17) 


.320545(10) 


-.938622(15) 


.29874(13) 


.19064(17) 


14 


-.962409(06) 


.320596(07) 


-.938641(07) 


.30017(11) 


.19183(09) 


18 


-.962434(14) 


.320660(09) 


-.938720(11) 


.29854(08) 


.19102(10) 


30 


-.962429(08) 


.320688(08) 


-.938717(05) 


.29677(09) 


.18925(11) 


oo 


-.96242365.. 


.32068921.. 


-.93872320 


.29677 


.18920 



Table 2c: Monte Carlo results for (3 = 0.50 
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